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.  *  ■ 

During  the  past  year  definitive  answers  were  achieved  for  several 
questions  raised  in  the  proposal  cf  February  1978.  These  include: 

1.  An  examination  of  the  relationship  between  the  scattering  matrix 
(the  Fourier  transform  of  the  scattering  operator)  and  the  integral  equations 
used  in  the  Singularity  Expansion  Method  (SEM)  established  that  only  the  com¬ 
plex  poles  off  the  axis  a^e  intrinsically  associated  with  the  scatterer.  While 
it  has  been  known  in  special  cases  that  those  on  the  axis  do  not  contribute  to 
the  field  (e.g.  Dolph,  "The  integral  equation  method  in  scattering  theory", 
Problems  in  Analysis,  Ed.  R.  C.  Gunning,  Princeton  University  Press,  1970. 

s 

An  AFOSR  Symposium)  ,v  this  appears  to  be  the  first  time  this  relationship  has 
been  clearly  exhibited.  Since  the  scattering  matrix  can  be  shown  to  be  analytic 
in  a  half-plane  containing  the  axis,  any  integral  equation  should  exhibit  the 
same  properties  for  this  region.  Those  of  SEM  fail  to  do  so  but  other  integral 
equations  which  do  can  be  given  in  all  cases. 

This  result  raises  the  interesting  and  as  yet  open  question  as  to 
whether  the  poles  of  the  integral  equation  of  SEM  in  addition  to  yielding  poles 
not  intrinsic  to  the  body  might  yield  only  a  subset  of  the  complex  poles  of  the 
scattering  matrix.  No  examples  where  this  occurs  are  yet  known  but  there  is  no 
proof  that  it  cannot  happen.  In  contrast,  if  an  integral  equation  using  the 
Green's  function  in  contrast  to  the  free  space  Green’s  function  is  used,  A.  G. 


1 


Tht  jt  fiQ  u/ar'lJ'X  fT/j?  cf  //j /  c /i  /pfe  f/?ccj  XfM{/  Co/Jfr/cx 
itigu  c*'  Expert 01*  ScQ/cr  q  //c{  /4 efer  SoQ^^/t 


S 1  xi Sf  ^  * '  f  *  t  \j -  ■  •  r  _  .  ~  -r  *  /  - 

Cl  PCCrS  /  /  CS  Q  t/Cf  C/tc  -fr<U"61tfeC  ^  f'Ctos' 


FINAL  REPORT  ON  AFOSR  GRANT  77-3358 


6/1/78  to  5/31/79 

*FO*ft»Tft‘  8  9-1  508 

Principal  Investigator:  Charles  L.  Dolph 


Table  of  Contents 


Item 


Page 


Charles  L.  Dolph,  Technical  Report  1 

Abstract  of  Tali:  at  American  Mathematical  Society  Meeting  5 

at  Biloxi,  Mississippi,  January  24-27,  1979 

On  the  Relationship  Between  the  Singularity  Expansion  6 

Method  and  the  Mathematical  Theory  of  Scattering 
by  C.  L.  Dolpii  and  S.  K.  Cho 

Review:  The  General  Theory  of  Eigenvibrations  in  the  39 

Theory  of  Diffraction,  by  N.  Voitovic,  B. 

Kacenelenbaum,  and  A.  Sivcv. 

Abstract:  A  Critique  cf  the  Singularity  Expansion  and  43 

Eigenmode  Expansion  Methods 

A  Critique  of  the  Singularity  Expansion  and  Eigenmode  44 

Expansion  Methods,  by  C.  L.  Dolph,  V.  Komkov  and 
R.  A.  Scott 

List  of  Publications  and  Curriculum  Vitae  for  A.  G.  Ramm  53 


C.  T.  Tai,  Zeros  of  the  Input  Impedance  Function  for  Thin 
Cylindrical  and  Prolate  Spheroidal  Antennas 


62 


Ramm  has  recently  shown  that  there  is  a  one-to-one  correspondence  between  the 
poles  of  this  Green's  function  and  the  zeros  of  the  eigenvalues  as  they  are 
used  in  the  Eigenmode  Expansion  theory.  This  will  be  elaborated  on  in  future 
work. 


2.  The  relationship  between  the  eigenvalues  of  the  integral  equations  of 
SEM  and  the  complex  eigenvalues  of  the  associated  partial  differential  equations 
—  whether  scalar  or  vector.  Tn  particular,  the  Integra1  orp'<»r-{r>r»s  of  SF.M- have 
at  most  two  eigenvalues  +  and  these  are  functions  of  the  at-most-denumerable 
number  of  complex  eigenvalues  of  the  associated  differential  equations. 

3.  These  results  were  first  presented  in  an  invited  address  at  a  special 
session  on  Integral  Equations  at  the  annual  meeting  of  the  American  Mathematical 
Society  held  in  Biloxi,  Mississippi,  January  24-28,  1979.  An  abstract  of  this 
talk  is  attached. 

Moreover,  because  of  the  unfamiliarity  of  many  SEM  workers  with  the 
mathematics  involved,  several  attempts  were  made  to  write  an  account  suitable 
for  the  IEEE  Transactions  on  Antennas  and  Propagation.  C.  E.  Baum  and  others 
felt  that  the  first  two  drafts  were  too  detailed  and  technical,  so  that  a  third 
attempt,  with  many  simplifications,  was  prepared.  The  resulting  caper,  joint 
with  S.  K.  Cho  of  the  Radiation  Laboratory  of  the  University  of  Michigan,  has 
been  submitted  to  the  above  journal.  It  is  entitled  "On  the  relationship  be¬ 
tween  the  singularity  expansion  method  and  the  mathematical  theory  of  scattering". 

While  a  copy  of  this  paper  in  its  submitted  form  is  attached  it  is  *ri 

doubtful  that  the  paper  will  appear  exactly  as  here.  Correspondence  with  C.  E. 
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Kentucky  seem  to  indicate  that  certain  sections  will  still  have  to  be  expanded. 
Knowledge  of  this  paper's  existence  seems  quite  wide  spread  in  view  of  the  num¬ 
ber  of  requests  that  have  been  received  for  preprints  from  various  SEM  workers. 

4.  As  indicat'  in  the  1979  proposal  it  was  felt  that  the  theory  of  non¬ 
self-adjoint  operators  in  Hilbert  space  "would  have  significant  implications  for 
that  part  of  the  electrical  engineering  formalism  known  as  the  Eigenmode  Expan¬ 
sion  Method  (EEM)  particularly  in  relevance  to  the  use  of  this  formalism  in  de¬ 
veloping  equivalent  circuits. 

A  survey  of  the  Russian  literature  revealed  that  a  significant  start 
in  this  direction  had  been  made  in  both  the  scalar  and  vector  cases  by  several 
Russian  workers.  In  particular,  papers  bv  A.  G.  Ramm  and  M.  S.  Agranovic  in 
the  scalar  case,  made  it  evident  a.  forteriori  that  the  EEM  formalism  based  on 
the  impedance  integral  equation  would  not  be  correct  in  general  and  thus  the 
work  on  equivalent  circuits  would  have  to  be  critically  examined. 

In  this  connection  C.  L.  Dolph,  with  considerable  help  from  V.  Komkov, 
reviewed  the  book:  The  Generalized  Methods  of  Eigenvib rat ions  in  the  Theory  of 
Diffraction  byVoitovic,  Kacenelenbaum  and  Sivov,  as  well  as  the  mathematical 
appendix  Spectral  Properties  of  Diffraction  Problems  by  Agranovic.  An  unedited 
copy  of  this  review  is  attached  but  it  is  expected  that  the  final  version  will 
appear  in  the  September  or  October  issue  of  Mathematical  Reviews  for  this  year. 

Because  of  its  importance  two  attempts  at  translation  are  underway. 

A  machine  translation,  arranged  by  Dr.  R.  Buchal  of  AFOSR,  and  an  American  Mathe¬ 
matical  Society  translation,  arranged  by  the  author  and  the  editors  of  Mathema¬ 
tical  Reviews. 
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As  the  attached  abstract  indicates,  much  of  the  above  material  was 
presented  in  a  talk  by  C.  L.  Dolph  at  the  International  Symposium  on  Recent 
Developments  in  Classical  Wave  Scattering  held  at  Ohio  State  University  June 
25-27,  1979. 

A  manuscript  by  C.  L.  Dolph,  V.  Komkov  and  R.  A.  Scott  entitled  "A 
Critique  of  the  Singularity  Expansion  and  Eigenmode  Expansion  Method"  will  ap¬ 
pear  in  book  form  in  the  conference  proceedings  to  be  published  by  Pergamon  Press. 

A  copy  of  this  manuscript  is  attached.  As  this  paper  contains  the  ^“vant  English 
translations  of  the  Russian  literature,  copies  of  its  bibliography  were  also 
distributed  (see  attached  copy)  .  Additional  copies  of  this  bibliography  will  be 
furnished  on  request. 

5.  With  the  aid  of  Dr.  R.  Buchal  of  AFOSR  and  Professor  F.  Gehring  and 
Lee  Zukowski  of  the  Mathematics  Department  of  the  University  of  Michigan  and 
many  others,  including  the  author,  A.  G.  Ramm,  formerly  of  the  University  of 
Leningrad,  was  successfully  brought  to  the  University  of  Michigan.  Since  his 
arrival  on  June  5th,  1979,  he  and  the  authors  have  been  actively  planning  and 
beginning  to  implement  further  research  in  the  areas  under  discussion.  A  de¬ 
tailed  plan  will  be  available  shortly.  A  list  of  publications  and  curriculum 
vitae  of  A.  G.  Ramm  is  attached. 
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Abstract  of  Talk  at  American  Mathematical  Society  Meeting 
Biloxi,  Mississippi,  January  Ik-21 ,  1979 


C.L.  Dolph, University  of  Michigan, Ann  Arbor, Michigan  48104«Predholc 
integral  equations, scattering  theory,  and  the  singularity  expanalor 
getbcd  (SUM). 

Since  197t  SSM  has  been  extensively  used  in  eleo tro magnetic  theory  a 
v  i  ident  from  the  review  articlo  by  C.li.  Baumt  Emerging  technology  for  tram 
d  nroad  band  analysis  and  synthesis  of  antennas  and  scatterera,"  Proe.  I. 
(y,  1976,  pp. 1 598-1616.  The  Fredholm  integral  equations  and  their  complex 
ngularities  used  in  SKM  will  be  discussed  and  i>£M  will  be  interpreted  in 
rns  of  mathematical  scattering  theory. 


ON  THE  RELATIONSHIP  BETWEEN  THE  SINGULARITY  EXPANSION 
METHOD  AND  THE  MATHEMATICAL  THEORY  OF  SCATTERING 


* 

C.  L.  Dolph 

Mathematics  Department 
The  University  of  Michigan 


S.  K.  Cho 

Radiation  Laboratory 
The  University  of  Michigan 


*Partially  supported  by  AFOSR-773358 
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1.  Introduction 


The  singularity  expansion  method  (SEM)  and  the  mathematical 
theory  of  scattering  have  been  extensively  developed  independently 
of  each  other  during  the  last  decade  or  so.  To  our  knowledge  the 
relationship  between  them  has  not  been  discussed.  In  this  paper 
the  implications  for  SEM  of  the  mathematical  scattering  theory 
will  be  explored.  As  a  result  some  important  conclusions  of  SEM 
will  be  seen  to  also  follow  from  scattering  theory  and,  at  the 
same  time,  certain  gaps  between  the  two  approaches  will  appear. 

The  SEM  was  introduced  by  C.  Baum  in  1971  [1]  as  a  technique 
for  solving  transient  electromagnetic  scattering  from  bodies  of 
finite  extent.  Since  then,  as  reference  to  [2],  [3]  and  [4] 

indicate,  not  only  has  it  been  widely  used,  but  an  extensive 
formalism  has  been  developed  for  it.  Succinctly,  the  SEM  is 
a  generalization  of  well-known  technique  of  linear  circuit  theory 
in  which  singularities  of  a  transfer  function  in  the  complex 
frequency  plane  are  used  to  determine  the  transient  response 
by  the  Heaviside  expansion  theorem.  More  details  of  the  pro¬ 
cedures  of  the  SEM  will  be  given  in  Section  2. 

Mathematical  scattering  theory  as  used  here  originates  in 
the  work  of  Lax  and  Phillips  [13,  14,  15]  for  the  scalar  case 
and  its  generalization  by  Schmidt  [23,  24]  for  the  electromagnetic 
case.  An  independent  development  of  the  theory  based  upon 
integral  equations  similar  to  those  employed  in  the  SEM  theory 
is  due  to  Shenk  and  Thoe  [25],  [26],  and  [27]  in  the  scalar  case 

and  to  Pyz'janov  [19],  [20],  and  [21]  in  the  electromagnetic  case. 


-  z- 


The  relevant  parts  of  this  theory  will  be  sketched  in  Section  3 
as  well  as  its  implications  for  SEM.  In  Section  4  additional 
aspects  of  the  SEM  formalism  will  be  commented  upon  and  some 
questions  raised  as  to  its  interpretation. 

2 .  The  SEM  Formalism  and  a  Mathematical  Interpretation 

In  view  of  several  descriptions  of  the  SEM  given  in  the 
reference  already  quoted,  only  a  very  brief  sketch  of  the 
relevant  formalism  will  be  given  here.  The  discussion  will 
further  be  limited  to  the  electromagnetic  case  for  convex  per¬ 
fectly  conducting  bodies. 

Maxwell's  equations  are  first  subjected  to  a  two-sided 
Laplace  transformation  with  respect  to  the  time  and  the  initial 
value  problem  replaced  by  a  time-independent  boundary  value 
problem  for  a  reduced  wave  equation  for  the  electric  or  magnetic 
field  subject  to  the  corresponding  boundary  condition  for  perfect 
conductivity  and  a  vector  radiation  condition.  A  representation 
for  the  scattered  field  is  then  sought  in  the  form  of  an  integral 
representation.  Application  of  the  boundary  condition  yields 
an  integral  equation  (or  in  some  instances  an  integral-differential 
equation)  whose  kernel  is  given  in  terms  of  some  free  space 
Green’s  function  or  dyadic.  The  unknown  in  the  integral  equation 
is  the  response  (often  a  current  density  on  the  scatterer)  induced 
by  the  incident  wave.  The  so-called  natural  frequencies 
{s^}  are  the  values  of  s  for  which  the  corresponding 
homogeneous  integral  equation  has  non-triviai  solutions.  At 


-  8  - 


-  J- 

such  a  value  of  s  ,  the  inverse  to  the  integral  equation  will 
not  exist  and  the  values  of  s  will  be  poles  of  the  inverse. 

The  res1:'  ting  residue  series  which  leads  to  coupling  coefficients 
will  not  be  discussed  here.  In  actual  practice  the  homogeneous 
integral  equation  is  replaced  by  a  matrix  equation  by  some  finite 
approximation  technique  such  as  that  of  the  method  of  moments. 

The  poles  are  then  sought  approximately  as  the  zeros  of  the 
determinant  of  this  matrix. 

The  formal  procedure,  leaving  aside  the  relation  betv;een  the 
solutions  of  the  integral  equation  and  the  matrix  equation,  has 
been  justified  mathematically  in  only  a  few  cases.  The  most 
complete  treatment  is  due  to  Marin  [17].  For  a  perfectly 
conducting  body  he  used  the  Franz  H-field  representation  and 
treated  the  resulting  magnetic  field  integral  equation  (MFIE ) 
for  the  induced  current  density;  namely, 

(I  -  ?L)-J  -  J  -  2  /  n  X  (VG  X  j)d  o  =-•  2(n  X  Hinc) , 

Here  n  is  the  outward  bound  unit  norma]  to  the  scatterer  and 
G  is  scalar  free  space  Green  function 

c-sLr“.r_'  I 

G  -  — - - - 

4tt  \r_~r_'  | 

Solutions  are  sought  in  the  Hilbert  space  of  elements  J  that 
are  tangent  to  T  and  square  integrable  on  I’  whore  T  is 
assumed  sufficiently  smooth.  opeei  f  ically  the  integral  equation 
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(I  -  4L2)-J  -  2(1  +  2L) * (n  *  HinC) 

is  shown  to  be  Fredholm  and  treated  by  the  method  of  Carleman. 

The  operators  (I  -  2L)  and  (I  +  2L)  are  shown  to  have 
inverses  for  the  same  values  of  s  and  consequently  these 
inverses  exist  whenever  this  last  integral  equation  has  an  inverse 
From  this  it  follows  that  (I  -  2L)  ^  has  a  pole  at  s  =  s 

n 

whenever  the  homogeneous  integral  equation  (I  -  2L) • j  =  0 
has  a  non-trivial  solution  at  s  -  s 

n 

The  poles  so  determined  lie  in  the  union  of  two  sets.  One 
set  consists  of  poles  corresponding  to  interior  resonances.  These 
arc  purely  imaginary  and  their  occurrence  is  due  to  the  method. 
They  do  not  contribute  to  the  scattered  field  and  hence  are  not 
intrinsic  to  the  scattering  problem.  The  other  set  consists  of 
poles  corresponding  to  exterior  resonances  and  they  arc  not 
purely  imaginary.  They  are  in  fact  eigenvalues  of  the  vector 
wave  equation 

2 

s 

V  x  V  x  H  +  — H  ~  0 
— n  2  _n 
c 

for  the  exterior  scattering  problem. 

The  Jlaim  that  these  poles  (exterior  resonances)  are  intrinsic 
the  scattering  problem  will  be  further  substantiated  by  their 
intorpreta tion  as  poles  of  the  scattering  matrix  in  Section  3. 
Marin's  conclusion  can  also  be  deduced  directly  from  the 

analytic  Fredholm  theorem.  In  its  general  form  the  theorem  is 
due  to  Steinberg  [28]  who  established  it  for  a  general  Banach 
space. 
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Recall  that  a  Banach  space  B  is  a  linear  space  in  which 
a  "distance''  is  defined  by  a  norm,  usually  denoted  by  [|  || 
and  it  is  complete  in  the  sense  that  every  Cauchy  sequence 
has  a  limit  (in  norm)  in  the  space.  In  general  there  is  no 
inner  product.  Let  ^f(B)  be  the  set  of  bounded  operators 
on  B  A  norm  in  ^(B)  can  be  defined  for  a  linear  operator 
A  by 


=  sup 
xGB 


11Ax!1 
II  x  II 


An  operator  A  in  ^(B)  is  called  compact  (or  completely 
continuous)  if  and  only  if  for  every  bounded  sequence  {x^} 
in  B  ,  {Tx^}  has  a  convergent  subsequence  in  B  . 

A  family  of  operators  r’(s)  depending  upon  s  is  called 
an  analytic  family  in  the  neighborhood  of  s^  if 


T  ( s )  =  l 


Vs's0> 


n 


v;here  convergence  ,s  in  the  operator  norm  and  where  is  in 

• 

The  theorem  can  be  stated  as  follows: 

If  T(s)  is  an  analytic  family  of  compact  operators  for 
s  G  K  ,  an  open  connected  subset  cf  the  complex  plane,  then 
either  (I  -  T(s))  is  nowhere  invertible  in  K  or  ((I  -  T(s)) 
is  meromorphic  in  K  .  At  a  pole,  the  equation  <J/  -  T(Sq)^  =  0 
has  a  non-zero  solution  in  B  . 
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If  B  is  a  separable  Hilbert  space  the  inner  product  can 
be  used  to  express  the  residue  in  terms  of  operators  of  finite 
rank.  That  is,  for  each  pole  there  will  exist  an  N  and  two 
sets  of  linearly  independent  functions 

such  that  the  corresponding  residue  can  be  expressed  as 

N 

l  (<fv*  * 

where  (f,g)  has  been  used  for  the  inner  product.  A  proof  of 
this  can  be  found  in  Reed  and  Simon  [22] . 

It  may  be  shown  directly/ or  Marin's  argument  can  be  used 
to  conclude, that  the  MFIE  is  a  Fredholm  integral  equation  of  the 
second  kind  having  an  integral  operator  which  is  compact.  It 
is  also  not  difficult  to  verify  that  the  MFIE  yields  a  family 
of  operators  analytic  in  s  .  Marin's  conclusion  then  follows 
from  the  Hilbert  space  analytic  Fredholm  theorem. 

However,  since  in  most  SEM  applications  it  is  necessary  to 
resort  to  numerical  methods,  it  is  mr re  appropriate  to  seek 
solutions  in  the  space  of  continuous  functions.  Such  a  space  here 
would  consist  of  the  Banach  space  of  all  continuous  J  tangent  to 
T  when  the  norm  is  taken  to  be 

||  j|j  =  max  | J(x)  |  • 
c  xST 

The  Fredhom  operator  of  the  MFIE  will  still  be  a  compact  analytic 
family  in  s  and  consequently  the  same  conclusions  about  the 
natural  frequencies  will  follow. 
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Unlike  the  papers  by  Marin  other  SEM  papers  are  difficult  to 
interpret  mathematically  since  neither  the  properties  of  the 
integral  equations  nor  the  spaces  in  which  solutions  are  sought 
are  specified.  Some  papers  also  employ  integral  equations  of  the 
first  kind  which,  it  is  claimed,  lead  to  approximate  solutions 
which  are  easier  to  handle  numerically.  It  should,  however, 
be  recognized  that  integral  equations  of  the  first  kind  cannot, 
unless  of  certain  special  types  like  Wiener-Hopf ,  be  used  to 
establish  existence  and  uniqueness  properties.  Integral  equations 
of  the  first  kind,  in  whatever  space  their  solutions  are  sought, 
are  not  a  well  set  problem  in  the  sense  of  Hadamard  so  that 
some  method  of  regularization  such  as  one  of  the  methods  of 
Tihonov  [30]  should  oe  used.  While  such  methods’ have  been 
used  in  other  electromagnetic  problems  [cf.  (8)]  they  do  not 
appear  to  have  occured  in  connection  with  SEM.  It  should  also 
be  emphasized  that  for  scatterers  of  finite  extent,  the 

Franz  E-field  representation  will  not 

lead  to  an  integral  equation,  no  matter  how  smooth  the  boundary 
fields  are  assumed  to  be. 

In  view  of  this  the  subsequent  interpretation  in  terms  of 
scattering  theory  will  be  limited  to  those  problems  involving 
Fredholm  equations  of  the  second  kind. 

3 .  The  scattering  operator  and  scattering  matrix. 

Several  considerations  have  influenced  our  decision  to 
restrict  further  discussion  to  scalar  rather  than  vector  wave 
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equations.  Not  only  are  all  the  essential  ideas  present -in  this 
simpler  case  but  all  proofs  of  the  assertions  we  will  make  are 
readily  available.  This  is,  unfortunately,  not  true  in  the 
other  case.  For  example  the  theorems  about  the  scattering  operator 
and  the  representation  of  the  scattering  matrix  are  stated  without 

s 

proof  by  Schmidt  [23,24].  He  indicates  that  they  follow  directly 
from  the  corresponding  results  contained  in  the  book  by  Lax-Phillips 
[13]  .  And  while  this  is  true,  we  do  not  expect  our  readers  to  have 
the  necessary  familiarity  with  the  Lax-Phillips  theory  to  find  this 
easy  to  do.  A  reference  to  Beale's  unpublished  thesis  [5]  ,  where 

many  of  these  details  are  carried  out,  will  soon  convince  one  of 
the  difficulties  involved.  Further,  as  noted  in  the  reviews  of  the 
papers  by  Pyz'janov  [19],  [20],  [21],  many  open  questions  concerning 

the  spaces  and  operators  used  by  this  author  are  also  not  specified 
and  thus  these  papers  are  open  to  the  same  criticism  as  most  of  the 
papers  in  SEM. 

Fortunately  a  straight  forward  extension  of  the  Shenk-Thoe 
theory  could  be  carried  out  for  s  ^  0  .  This  would  involve 
replacing  the  representation  for  the  scalar  field  developed  by 
Werner  [31]  and  used  by  them  by  Werner's  subsequent  representation 
of  scattered  electric  field  as  given  in  [32] .  The  Shenk-Thoe 
procedure  could  then  be  followed  step-by-step  and  the  con¬ 
clusions  reached  by  Schmidt  obtained  in  this  way  just  as  the 
conclusions  of  the  Lax-Phillips  theory  can  be  reached  by  the 
Shenk-Thoe  method  in  the  scalar  case. 

While  Laplace  transformations  in  time  have  so  far  been  used, 
henceforth  they  will  be  replaced  by  an  equivalent  Fourier  transform. 


The  time  convention  used  will  be  that  of  most  papers  in  physics 
and  under  it,  the  upper  half  plane  replaces  the  right-hand  half 
plane  as  the  domain  of  analyticity.  This  is  also  the  convention 
used  by  Shenk-Thoe  and  its  use  allows  one  to  state  spectral 
properties  in  the  usual  form.  It  amounts  to  introducing  k  =  is  . 
It  is  not  that  of  Lax-Phillips  and  Schmidt  where  "i"  would 
be  "-i"  .  Thus  outgoing  waves  here  would  be  incoming- 
in  their  sense  and  vice-versa. 

Let  ft  be  an  unbounded  free  space  (ft  C  En ,  n  =  2,3) 
exterior  to  a  sufficiently  smooth  convex  body  with  the  surface  T  . 
The  time-dependent  acoustic  scattering  problem  seeks  the  solution 
for  the  wave  equation  for  the  scattered  field  us 

1  g2 

(A  -  ~2  — 2^  uS(r't)  =  0  ,  r  G  f!  ,  t  >  0 

c  3 1 

subject  to  the  initial  conditions 

us  (r , 0)  =  f 1 (r)  , 

t  =  0  ,  r  6  ft  (3.1) 

~  us  (r  ,  0)  =  f  2  (r)  . 

and  one  of  the  following  three  boundary  conditions  for 
t  >  0: 


Dirichlet : 

us(rQ,t)  =  g1(r0»t)  , 

Neumann : 

3n0  U  (r0't)  =  92(r0't)' 

r  p  r 
0  1 

Robin : 

{3n0  +  a)  u  “  ?3(r0 

,  t)  , 

-  15  - 
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where  g-pj—  =  n^-V^,  ng  denoting  the  unit  normal  vector  at 
the  observation  point  on  T  directed  into  ft  . 

In  the  scattering  theory,  one  views  these  problems  as 
perturbations  of  the  free  space,  or  unperturbed  problem  (with¬ 
out  the  scattering  obstacle)  such  that,  for  large  time  |t(, 
the  perturbed  (or  scattered)  solutions  become  asymptotically 
equal  to  the  free  space  solutions.  The  scattering  operators 
which  are  thought  to  contain  all  the  observable  information  about 
the  nature  of  the  scatterers  represent  measures  of  these 
perturbations.  Let  U(t)  denote  operators  which  relate  the 
states  at  t  =  0  to  those  at  t  >  0  in  the  scattering  problems. 
Similarly,  let  UQ(t)  denote  the  corresponding  operators  in  the 
free  space  problems.  If  f  =  (f^f^  denotes  the  initial  data, 
then  U{t)f  solves  a  scattering  problem  in  the  sense  that  it 
represents  the  state  at  time  t  >  0  ,  evolved  from  the  initial 
state.  Since  the  scattered  solution  is  assumed  to  be  asymp¬ 
totically  equal  in  large  |t|  to  the  free  space  solution,  one 
expects  that  there  exist  initial  states,  fQ  +  ,  in  the  fre .  space 
unperturbed  problem  such  that  Ug (t)  f^  +  are  asymptotically 
indistinguishable  from  U(t)f  as  t  •>  +  «*.  This  is  made  precise 
by  introducing  an  appropriate  Hilbert  space  with  the  energy  norm 

Ilf!!  =  /  MfJ2  +  I  f  2 1  2J  dv 

ft  x  * 

where  f  ^  and  f^  are  compactly  supported.  For  this  space, 
the  operators  U(t)#Ug(t)  are  one-parameter  groups  and  unitary. 

The  asymptotic  behaviors  stated  above  are  now  expressed  as 
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II U  ( t )  f  -  UQ(t)  f  Q—  H  -»■  0  ,  as  t  ->  -00  , 

||u(t)f  -  U0(t)f0+  |!  ->0  ,  as  t  +00  . 

Since  U(t),  U^t)  are  unitary  and  form  a  group,  the  above 
expressions  may  be  written  as 

||u(~t)  'U(t)f  -  UQ(t)  fQ+]  ||  0,  as  t  -*■  +-  ; 

or 

||  f  -  U(-t)U0{t)f0+||  -*•  0  ,  as  t  -  +»  . 

The  operators 

W+(t)  =  lim  U(-t)U0(t)  , 

—  t  ->  +°°  u 

exist,  and  are  called  the  wave  operators  as  is  well  known  in 
the  physics.  V.'  (t)  ,  for  instance,  compares  the  asymptotic 

free  space  solution  UQ(t)fQ-  with  the  asymptotic  scattered 
solution  U(t)  fQ+  .  See  the  Figure  1  below. 


The  map 


I - > 


V 


(3.2) 
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is  called  the  scattering  operator.  That  such  operators  exist 

and  are  unitary  in  the  Hilbert  space  can  be  established.  As 
evident  from  Figure  1/ 

i  =  w+~l  w.  • 

Let  denote  the  n-dimensional  (n  =  2,3)  Fourier 

transformation  defined  by 

tJ-(f)  =  f  (k)  =  (27r)"n/2  /  f(r)  e"*1^ ' r  d’V  .  (3.3) 

—oo 

The  scattering  matrix  is  then  defined  by 

•  0.4) 


Next,  consider  an  exterior  homogeneous  boundary  value 
problem 


(A  +  k2)  d+(r)  =  0  ,  r  6  0  , 
|=0,  F0  6  F  ' 


(3.5) 


Here  y  is  written  for  the  Dirichlet,  Neumann  or  Robin  boundary 
condition. 

Then , 

<fr+(r)  =  eX*'r  +  v+  (r  ,k)  (3.6) 

are  evidently  the  solutions  of  the  problem.  Although  they  are 
not  in  1, 2  (ft)  »  <j)+(r)  arG  improper  eigensolut  ions .  v+(r,k) 

are  the  solutions  of  the  boundary  value  problem 
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(A  +  k2)  v+(r,k)  =  0,  r  6  . 


Y[v+  (rQ,k)  + 


ik-  r 
e 


(3.7) 


v+  (v_)  is  called  an  outgoing  (incoming)  diffracted  plane 
wave;  <{>+(r)  are  called  distorted  plane  waves. 

If  a  generalized  transformation  is  defined  by 


-%  +  (f)  =  lim  {2TT)“n/2  /  f(r)  <f>+  (r  ,k)  *dnr  , 

—  Q 

P 


(3.8) 


where  Qp  =  { r  S  Q  :  Jr|  <  p}  for  some  sufficiently  large  P 
and  where  *  denotes  complex  conjugate,  then  the  scattering  matrix 
S  takes  a  simple  form 

s  =  %  ST1  , 

~r  — 

which  is  equivalent  to 

f  =  f  . 

Th  the  scattering  matrix  maps  <p_  (r,k)  into  <f>+ ( r,kj  .  This 

form  is  particularly  convenient  for  time-independent  scattering 
theory.  Since  k  will  be  continued  to  the  complex  plane,  it 
is  necessary  to  replace  the  usual  radiation  condition  by  a 
representation  based  on  Green's  formula.  For  v+(r,k)  this 
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takes  the  following  form:  for  outgoing  waves  associated  with 
the  plus  sign  and  incoming  waves  associated  with  the  minus  sign, 


v±  (r,k)  =  /  [v±  (.r '  ,k)  g|r  G1  (kR)  -  C1  (kR)  ~  v± 

(1) 

where  G”  (kR)  =  ±  |  (^r)  (n  2)/2  H  (n-2)  /2  {kR)  '  n 

R  =  I?-?' I ,  r  e  a  ,  F*  e  r  . 


(r',k)]dA*  , 


=  2,3 


(3 .  S) 


This  representation  implies  that  the  corresponding  integral 
of  the  right-hand  side  taken  over  a  large  sphere  tends  to  zero 
as  the  radius  of  this  sphere  tends  to  infinity.  It  can  be  shown 
that  this  in  turn  is  equivalent  to  the  usual  radiation  conditions 
for  k  ^  0  and  0  <  arg  k  _<  tt  . 

For  k  real,  r  =  |rj  ,  0  =  r/r  ,  w  =  k/k,  k  =  |k|  ,  and  for 

| r  j  »  [ r ’  !  on  obtains  the  asymptotic  form 

_  _  _  _  ±ikr 

v+(r,k)  =  v+(r0,kui)  =  -  --  j^iyy/2  fs+(’9'k/w)  +  0(1)]  .  (3.10) 


The  functions  s+  in  this  last  formula  are  called  radiation 
patterns  or  transmission  coefficients.  In  terms  of  them  the 
scattering  matrix  takes  the  form 


S  (k)  h  (8)  =  h  ( 9  )  + 


( ik/ 2tt )  (n  1)/2  / 


OJ 


h (w) s_ (-0 ,k ,w) *dS 
-  1  w 


(3.11) 


for  h  ( 0 )  e  ( Sn  1)  ,  Sn  1  being  a  surface  of  the  unit  ball 

in  n-dirnensional  space.  The  above  integral  operator  is 
compact  for  k  >  0  . 
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This  representation  may  be  derived  from  the  scattering  operator 
[26]  or  it  may  be  verified  by  first  postulating  it  and  then 
invoking  some  form  of  the  radiation  condition  and  Hellich's 
uniqueness  theorem,  [13,  15]. 

The  representation  formulas  for  outgoing  (incoming)  waves 
v+  imply  that 


s+  ( 8  /  k ,  cj) 


=  ±i 


8  (k) 

2tt 


/  fv,  (y,k,w) 

r  - 


.  .9  +ik  y  *e 

3n 1 


where  8 (k) 


+  ik  v  *  8  3  -  _ 

e  3n-rv±(y,l:,u)]day 

_  ..  (^) 

(+  — )  .  n  =  2 , 3 ,  k  >  0  . 


(3.12) 


This  relation  establishes  a  one-to-one  map  between  the 
analvtical  oroperities  of  s ,  and  of  v,  .  Since  it  can  be 
shown  that  s_  (-7)  , k  ,  w)  *  =  s^_  (-9  ,  k*  , -«)  it  follows  that  if  the 

T 

solution  v+  can  be  constructed  and  shown  to  have  a  meromorphic 
continuation  from  the  positive  real  axis  to  the  entire  complex 
plane  if  n  is  odd  and  to  the  entire  logarithmic  Riemann 
surface  { k  /  0  ,  <  arg  k  <  °°  if  n  is  even  ,  it  will  follow 

that  the  same  is  true  for  the  scattering  matrix.  It  also  follows 
from  the  Shenk-Thoe  theory  that  the  scattering  matrix  is  unitary 
for  each  positive  k  . 

Shenk  and  Thoe  construct  the  solution  for  v,  for  a  more 
general  problem  in  which  surface  T  may  have  several  components 
and  admit,  a  perhaps  empty,  decomposi tion  as  the  union  of 
ri  and  r 2  -  On  Robin  boundary  conditions  are  imposed 

(a  =  0  is  allowed)  while  Dirichlet  boundary  conditions  are 
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imposed  on  r .  The  scattered  field  is  sought,  following 
Werner  [31]  ,  as  a  superposition  of  a  single  layer,  a  double 
layer  and  a  volume  potential  and  the  resulting  system  of  Fredholm 
integral  equations  of  the  second  kind  is  shown  to  be  uniquely 
solvable  for  Im  k  >  0  in  an  appropriately  constructed  Banach 
space.  Steinberg's  theorem  is  then  invoked  ^to  establish  the 
analytic  properties  of  the  scattering  matrix.  This  method, 
while  quite  general,  would  not  be  useful  for  numerical  computation 
because  of  the  presence  of  the  volume  potential  and,  to  our 
knowledge,  it  has  never  been  employed  for  this  purpose. 

In  contrast, for  the  Dirichlet  and  Neumann  problem  alternate 
representations  using  a  complex  combination  of  single  and  double 
layer  potentials  not  only  can  be  used  to  establish  a  1-1  map  , 
but  is  more  useful  for  numerical  calculations. 

The  represer.tat  i  ons  are  due  to  Brakhage  and  Werner  [7]  and 
to  Lois  [1^]  and  take  the  form: 


v+(r,k)  =  /  [—■  --  -t  ixG+(kR)j  ||)(y)da  , 

r  y  y 


v+ (r ,k) 


4 

/  [G+(kR)  +  ix  ]  ifi(y)dax 


r 


(3. 

(3- 


Here,  i  ~  1  for  Re  k  >.  0  and  x  =  -1  for  Re  k  <  0  .  The 
resulting  integral,  equations  for  a  plane  wave  incidence  arc 

<Mx)  +  2  /  +  i T G+ )  (y )  da  =  -2e  x 

f  6  y  y 


(3. 


a  nd 


«M"z)  -  2  /  [' 


3_G_ 

In 


IT 


X 


3  n  3  n 
x  y 


<My)do  - 


3  n 

x 


-ik  •  x 


(3 


13a) 


13b) 


14a) 


14b) 


2  3  - 


These  equations  have  analytic  compact  integral  operators 
in  the  Banach  space  of  continuous  functions  and  possess  a  unique 
solution  for  Im  k  >  0  and  are  meromorphic  elsewhere.  They 
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therefore  also  provide  a  1-1  map  betv;een  the  transmission 
coefficient,  and  the  solutions  to  the  reduced  wave  equation. 

It  should  be  noted,  however,  that  because  of  the  higher  order 
singularity  in  the  Neumann  problem  it  is  necessary  to  use  a 
method  of  regularization  to  obtain  these  results.  This  is 

carried  out,  follow; ng  the  method  of  Leis  [16]  ,  in  detail  in 

/ 

Kussmaul  [12]  where  the  cylinder  problem  is  treated  numerically. 
The  corresponding  numerical  treatment  for  the  cylinder  Dirichlet 
problem  is  given  by  Greenspan  and  Werner  [10] . 

If  SEM  were  specialized  to  this  case  the  corresponding 
assumptions  would  be  that 


v 


+ 


G+6 (y ) doy 


for  the  Dirichlet  problem  and 


(3.15) 


v 


/  G+t(y)day 


(3.16) 


for  the  Neumann  problem.  The  corresponding  integral 
equation  weald  be,  respectively, 

<5  (x)  +2  /  6  (y) da  -  -2  (3.17) 

T  y  y 

y(x)  -  2  J  G+y(y)da  =  2  e~lE'X  (3.18) 

r  x  x 

Since  as  shown  in  Kupradze  [11]  these  representations  do  not 
provide  a  unique  solution  to  the  corresponding  boundary  value 
problem  for  the  reduced  wave  equation,  the  non-trivial  solutions 
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at  (kn)  of  the  homogeneous  equation  are  not  in  1-1  correspondence 
with  the  poles  of  scattering  matrix.  The  scattering  matrix  has 
no  poles  in  the  upper  half-plane  including  the  real  axis  and 
thus  does  not  contain  the  corresponding  to  interior 

resonances.  Nor  do  these  occur  in  the  expression  for  the 
scattered  field.  A  direct  proof  of  this  last  fact  for  the 
Dirichlet  problem  can  be  found  in  Dolph  [9] .  The  complex  » 

as  may  be  seen  by  direct  computation,  to  correspond  to  outgoing 
solutions  of  the  reduced  wave  equation  are  complex  eigen¬ 
values  of  the  exterior  scattering  problem  and  consequently 
correspond  to  poles  of  the  s-matrix.  These  are  method  independent 
and  intrinsic  to  the  scattering  body  in  contrast  to  the  interior 
resonances  which  arise  only  because  of  the  method. 

Because  these  presentations  fail  to  furnish  a  1-1  map,  it  is 
also  conceivable  that  there  may  be  other  poles  of  the  scattering 
matrix  not  given  by  solutions  of  the  homogeneous  integral  equations 
used  in  SEM.  There  are  no  known  examples  in  which  this  occurs, 
but  the  problem  must  be  considered  open.. 

More  explicitly  the  complex  poles  of  the  scattering  matrix 
are  in  1-1  correspondence  with  the  solutions  of  the  boundary 
value  problem 


(A  +  kn^Wn^r'^n^  =  r  6  ^ 

Ywn(70'kn}  “  0  '  r0  6  r  * 
The  complex  poles  kfi,  n  =  1,2,3,...  , 


(3.19) 


26  - 


kn  =  <k0*n  -1  k;  '  kn  >  0  ' 


0  >  Im  (k^)  >  ImCk^)  >  ImCk^)  > 


lead  to  the  representation  of  the  solution  of  (3.1)  in  the 
form 

_  _  “  -ik  t  __ 

u  (r,t)  =  1  a  e  w  (r) .  (3.20) 

n  n 

n=l 


While  the  functions  w  (r)  grow  exponentially,  in  odd  dimensions 
the  causality  implies  that  solution  are  zero  for  time  t  <  ~  >- 
c  being  the  speed  of  the  wave  propagation;  after  this  time, 
the  solution  decays  in  time. 

Lax-Phillips  [13]  proved  that,  under  the  Dirichlet  boundary 
condition  for  n  =  3  ,  the  scattering  operator  uniquely  determines 
the  scatterer.  There  are  few  theorems  about  the  locations  of  the 
complex  poles  of  the  scattering  except  for  those  located  on  the 
negative  imaginary  axis  (in  the  SEM,  on  the  negative  real  axis) 
of  the  complex  wave-number  plane.  Thus,  Lax-Phillips  showed 
that,  under  Dirichlet,  or  Neumann  boundary  condition  in  3-dimension, 
if  a  =  Im(k)  ,  Re(k)  =  0  and  if  the  scatterer  contains  a  sphere 
of  radius  ,  and  is  contained  in  a  sphere  of  radius  R2  , 

then  the  number  of  poles  N (o )  is 


lim  inf 

0-)-co 


N  (o) 
_  2 


>  -(^)2 
-  2  Y  0 


t 


v/hile,  if  the  scatterer  is  star-shaped, 
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lim  sup 
a-»“ 


N(o) 


1  R 

I(_2 

2  Y0 


I 


Y 0  =  0.66272  . 

J.  Beale  [6]  extended  the  above  results  to  all  n  ,  i.e., 
n  =  2,3  ,  including  the  case  of  Robin  condition.  In  particular,  for  i 
he  showed  that  the  scattering  matrix  has  at  most  a  finite  number 
of  poles  on  the  purely  imaginary  (purely  real)  axis  when  a 
Fourier  (Laplace) transformation  is  used.  Furthermore,  for  n  =  2  , 
if  the  boundary  condition  is  either  Dirichlet,  or  Neumann,  then 
there  are  no  poles  on  the  axis.  The  poles  on  the  axis  are,  of 
course,  purely  decaying  modes.  The  complex  poles  of  the  scattering 
matrix  which  are  also  the  complex  eigenvalues  of  in  an 

exterior  homogeneous  boundary  value  problem  for  the  Helmholtz 
equation  are  associated  with  the  physical  nature  of  the  scattering. 

Finally,  ve  will  quote  an  example  due  to  Shenk-Thoe  [27] 
to  show  that  the  complex  poles  of  the  scattering  matrix  are 
precisely  those  exhibited  in  the  separation  of  variable  solution. 

Let  Sn  ^  ,  n  =  2,3  ,  denote  the  unit  ball  in  En  .  Then, 
decompose  L2(Sn  ^)  into  finite  dimensional  subspaces  of 

spherical  harmonic  of  degree  l  .  Each  of  these  spaces  is  an 
eigenspace  of  the  operator  S(k),  the  scattering  matrix,  and  its 
eigenvalue  is 
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_  JB+i 
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where  a  denotes  the  radius  of  sphere,  or  circle.  The  .  . 

fact  that  the  numerator  and  the  denominator  of  the  eigenvalue 

of  S(k)  above  are  complex  conjugate  of  each  other  exhibits 

the  unitarity  of  S (k)  for  k  >  0  .  The  poles  of  S (k) 

are  the  roots  of  (ka)  which  occur  when  k  is  complex  . 

p+x, 

In  particular,  for  n  =  2  ,  consider  an  arbitrary  continuous 
function  on  the  circle  of  radius  a  .  Expanding  it  in  Fourier 
series  ,  one  obtains 


S(k)  (l 
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(3.21) 


where  (a,  <J>q)  denotes  a  point  on  the  circle.  If  k  is 
complex  ,  S(k)  has  complex  poles  in  the  lower  half-plane 
of  the  complex  wavenumber  plane. 

Since  the  complex  singularities  determined  by  integral 
equations  depend  only  upon  the  homogeneous  integral  equation, 
it  suffices  to  consider  the  simple  boundary  value  problem 


(A  +  k2)uS{r)  =  0  ,  r  6  ft  . 


(3.22) 
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This  has  the  known  separation  of  variable  solution  given  by 


<»  J  (ka)  , ,  .  .  , 

u5(F)  =  -  l  <-i)n  -  Ii'1)(kr)elr",>  .  (3.23) 
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This  solution  exhibits  the  same  complex  singularities  as  the 
scattering  matrix.  The  Brakhage-Werner  representation  f7] 

uS(?)  =  f2J  iT)  H^1}  (kR)u  .  (3.24) 

leads  to  the  integral  equation: 

V  (a,4>0)  +  — ^  /  (aTklT  ~  12t>  (kR0)y  d$  '  (3.25) 

-ika  cosij). 

=  -2  e  u 

The  unique  solution  of  this  integral  equation  is 

in<f,0 
e 

-  .  (3.26) 

-  it Jn (ka) ] 

Again  the  complex  singularities  are  the  same  as  those  of 
the  scattering  matrix. 

The  SEM  type  of  representation  for  this  problem  would 
be 

uS(7)=  J  Fife')  Ho1}  (kR)  v  (a, <|,')  dem¬ 

and  the  corresponding  integral  equation  would  be 

v(a>V  +  t/,?®  Ho1>(kRo,v  d»' 

-ika  cos  <j>n 
=  -2  e  U 

with  a  solution  given  by 
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v  (a ,  4>  0 ) 
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irka 


l 


(-i)n  1Jn(ka)  in4>0 
J^(ka)H^1)  (ka) 


(3.29) 


Here,  in  addition  to  the  complex  singularities  given  by  the 

scattering  matrix,  the  zeros  of  (ka)  occur  at 

those  values  of  k  which  correspond  to  the  eigenvalues  of  the 


corresponding  interior  homogeneous  Neumann  problem.  They  are  an  exampl 
of  the  interior  resonances  of  SEM  and  they  do  not  appear  in  the 

separated  solution  for  the  scattered  wave.  They  are  clearly 
seen  to  be  method  dependent. 


4 .  Additional  Comments  on  the  SEM 

Eigenmode  expansions  and  other  spectral  properties  have  been 
employed  in  the  development  of  the  SEM  formalism  by  C.  Baum 
[2,3,4]  and  L.  Marin  and  R.  Latham  [18].  While  these  have  been 
used  only  formally,  even  the  formal  expressions  are  often  meaning¬ 
less.  For  example,  the  left  and  right  eigenvalues  are  defined 
formally  by 

<  r(r,f';s);Rs(r';s)  >  =  X(s)Rg(r;s)  , 

<  Lg  (r;s)  ;_T(r,r ,  '  ;s)  >  =  X(s)Lg(r;s) 

where  _T(r,r';s)  /represents  some  intec  ral  operator  involving 
a  free  space  dydadic  kernel.  Then,  the  integral  equation 
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<Cr(r,r';s);  J  (r '  ;  s)>  =  I  (r ;  s) 
has  a  formal  solution  given  by 

i  <L  (r;s);  I(r';s)  > 

J(r;s)  =  l  — -H__ - - -  r  (r;s)  . 

p  Xp(.$)  <Lp(r;s);  Rp(r*;s)>  p 

Since  neither  properties  of  the  integral  equation,  nor 
the  space  in  which  its  solutions  are  sought  are  specified, 
it  is  difficult  to  comment  on  the  validity  of  this  procedure. 

If  the  space  were  a  Hilbert  space  of  square  integrable  functions 
and  the  integral  equation  were  of  the  first  kind  with  a 

hermitian  symmetric  Hilbert-Schmidt  kernel  then 
under  a  strong  convergence  condition  it  is  valid.  (Cf.e.g., 
the  discussion  in  Tricomi  [29]). 

An  alternate  and  more  promising,  interpretation  would  be 
in  terms  of  non-self  adjoint  operators  as  they  are,  for  example, 
discussed  in  Ramm  [21.1]  .  This  interpretation,  however, 

requires  a  knowledge  of  dissipative  operators  and  root  spaces 
and  consequently  will  not  be  further  commented  upon  here. 

It  appears  that  much  of  the  doubtful  SEM  formalism  arises 
whenever  a  formal  limit  of  a  matrix  relation  is  taken.  While 
all  matrices  have  eigenvalues  and  associated  invariant  spaces, 
this  is  not  true  of  their  limit,  integral  operators,  and,  even 
if  the  integral  operator  does  have  a  eigenvalues  there  may  not 
.  be  enough  to  span  the  range  of  the  operator.  As  a  case  in 
point,  the  homogeneous  MFIE  has  only  one  eigenvalue,  namely  "-l" 
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but  this  value  can  arise  for  many  different  values  of  k-those 
corresponding  to  both  interior  and  exterior  resonances.  While 
these  k's  corresponding  to  exterior  resonanaces  are  eigenvalues 
of  the  exterior  scattering  problem  for  the  vector  wave  equation, 
they  are  not  eigenvalues  of  the  MFIE. 

Since  there  is  a  rapidly  developing  theory  of  equivalent 
circuits  based  on  the  SEM  formalism,  this  formalism  should  be 
carefully  evaluated  and  only  that  part  of  it  for  which  a  mathe¬ 
matical  basis  can  be  furnished  should  be  used.  We  hope  to  return 
to  this  subject  in  the  future. 
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REVIEW 


Unedited  draft  -  The  general  method  of  eigenvibrations 

in  the  theory  of  diffraction. (by  N.  Voitovic,  B.  Kacenelenbaum, A. Sivov 

Spectral  prop,  of  diffractions  problems  (by  M.S.  Agranovic  ) 

NAUK  Moscow,  1977. 

The  problems  of  diffraction  and  scattering  of  electro¬ 
magnetic  or  ccoustic  waves  have  been  treated  at  different  levels 
of  physical  evidence  or  intuition  and  of  mathematical  rigor.  In 
some  engineering  publications  "the  analysis"  consisted  of  little 
more  than  formal  manipulation  without  much  worry  about  the 
spaces  in  which  the  operators  act,  about  completeness  or  con¬ 
vergence.  At  the  other  end  of  the  scale  sophisticated  mathematical 
theories  have  been  developed  which  so  far  failed  to  permeate  to 
the  level  of  physical  applications.  The  monograph  under  review 
is  a  splendid  example  of  the  intermediate  "interdisciplinary 
approach,  attempting  to  open  up  the  communication  between  the 
mathematicians  and  practical  engineers  working  on  problems  of 
communication  theory  In  the  opinion  of  the  reviewer 

J 

their  approach  appears  to  be  the  most  promising  way  of  overcoming 
the  vocabulary  and  traditional  formalistic  differences  in  approach 
to  the  subject.  The  first  part  outlines  some  eigenvalue-eigen¬ 
function  problems  of  electromagnetic  theory,  particularly  the 
£ -technique  which  has  been  pioneered  by  Voitovic,  Kacenelenbaum, 
and  Sivov  The  "generalized  method"  consists  of  representing  the 
solutions  in  terms  of  eigenfunctions  of  a  homogeneous  problem  in 
which  the  eigenvalue  is  not  the  frequency  (or  a  root  of  the  fre¬ 
quency)  but  some  other  parameter  of  the  system,  such  as  for  example, 
the  dielectric  constant  of  a  reference  system  occupying  the  same 
volume.  It  is  not  necessary  for  each  term  of  the  series  to 
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satisfy  the  boundary  conditions  (as  is  the  usual  practice  in 
applying  the  Rayleigh-Ritz  or  Galerkin  procedures) .  Only  the 
whole  series  satisfies  the  imposed  boundary  conditions.  Part  two 
considers  the  significance  of  the  boundary  conditions.  The 
conditions  of  impedence  type  (the  w  techniques) ,  the  adjoint 
boundary  conditions  (the  p-techniques) ,  conditions  given  at 
infinity,  the  metallic  surface  conditions  are  some  of  the  topics 
discussed.  Chapter  3  covers  the  variational  techniques.  Using 
the  e-approach  they  considered  the  system 

k*eu  =a  0  u  y  + 

r  ?  * 

*“=*0  B  y~ 


u+ 


du+ 


du~  J 

**  Is 


“Is  — 0. 


Considering  e  as  a  (fixed)  constant,  the  authors  look  for  the 
stationary  behavior  of  the  following  functional 
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or  of 

L(u)  /  (^u) 2  dV  -  k2  /_  u2  dV  -  k2s  f+  u2  dV. 
v  v  v 

Alternately  setting  L(u)  =  0.  They  consider  the  functional 

E(u)  of  the  e-technique 

\  (Vi«)*  dV  -  k1  \  u*dV 

- r - ~ - • 

k*  \  u'dV 

Y+ 


Variants  of  these  techniques  are  applied  to  a  variety  of  problems 
in  Chapter  4.  Metallic  tuning  devices,  wave  guides  constitute 
typical  examples  of  applications.  A  separate  chaptep, 

^pectral  properties  of  diffraction  problems  (123  pages  plus  biblio- 
graphy)  ,  v/rirten  by  Agranovic  is  modestly  called  an  appendix.  It 
fully  deserves  a  separate  review.  It  contains  the  introductory 
definitions  concerning  non-selfadjoint  operators  in  Hilbert  spaces, 
basis  and  biorthogonal  systems,  Sobolev  spaces,  and  elliptic 
operators.  He  discusses  some  properties  of  elliptic  pseudodiffer¬ 
ential  operators  and  of  elliptic  boundary  value  problems.  This 
preliminary  discussion  leads  to  theorems  concerning  the  ellipticity 
of  the  boundary  value  problems  in  diffraction  theory,  and  the  order 
of  the  corresponding  operators  and  of  their  selfadjointness. 
Problems  in  diffraction  theory  are  reduced  to  the  study  of  spectral 
properties  of  certain  pseudodifferential  operators  theorems  on 
disipative  property  of  such  operators  are  proved.  The  author’s 
s-technique  is  rigorously  restated.  Invertibility  (which  was 
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required  in  certain  steps  of  Chapter  2  is  proved  as  a  theorem. 
Some  properties  of  the  s-technique  are  closely  related  to  the 
general  techniques  initiated  by  M.G.  Krein.  The  appendix  does 
not  discuss  the  details  of  the  variational  approach  of  the 
authors.  Further  research  is  needed  to  resolve  many  outstanding 
issues  raised  by  the  authors  approach  to  diffraction  theory. 
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A  CRITIQUE  OF  THE  SINGULARITY 
EXPANSION  AND  EIGENMODE  EXPANSION  METHODS 


C.  L.  Dolph,  V.  Komkov  and  R.  A.  Scott 
University  of  Michigan,  Ann  Arbor,  Michigan  48109 


SUMMARY 

The  authors  outline  some  available  theoretical  techniques 
interpreting  the  complex  singularities  of  the  S  matrices,  and 
po-nt  out  various  difficulties  which  so  far  have  prevented  a 
formulation  of  a  consistent  and  vigorous  theory.  Some  directions 
for  future  research  are  suggested.  These  are  primarily  concerned 
with  the  singularity  expansion  and  the  eigenmode  expansion  methods. 
In  particular,  some  recent  results  of  the  Russian  school  pertaining 
to  the  theory  of  non-self-adjoint  operators  are  summarized. 
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The  role  of  complex  singularities  in  scattering  problems  is  a  rich  and  diverse  one; 
[see  Dolph  and  Scott  (1)].  Baum  (2)  has  been  of  the  prime  movers  in  the  United  ' 
States  of  the  subjects  of  the  title.  He  has  been  helped  in  his  efforts  by  the  work 
of  Marin  (3)  and  Tesche  (4) .  In  the  Soviet  Union  their  counterparts  are  Voitovic,  ' 
Kacenelenbaum  and  Sivov  (5).  There  appears,  however,  to  be  an  important  .difference  j 
between  the  two  countries  in  that  in  the  USSR  mathematicians,  beginning  with  A.  G.  ? 
Raima,  carefully  examined  the  formalism  of  the  approaches.  Such  efforts  have  j 

culminated  in  a  book  with  an  interdisciplinary  spirit  entitled:  The  Generalized 
Methods  of  Eigen  Vibrations  in  the  Theory  of^  Diffraction  (5).  The  book,  which  also, 
contains- an  extensive  appendix  by  Agranovic- entitled :  Spectral  Properties  of  ; 

Diffraction  Problems,  has  been  reviewed  by  the  first  author  (with  considerable  help 
from  the  second  author) .  This  review  should  appear  shortly  in  Mathematical  Uevievrs. 
and  an  effort  is  underway  to  have  the  text  itself  translated  into  English.  It  is  ; 
fortunate  that  much  of  the  pertinent  Soviet  work  has  appeared  in  English  translaticn 
in  the  Journal  of  Radio  Engineering  and  Electron  Physics.  The  readers  'attention  j 
should  be  called  to  the  two  papers  of  Ramm  [(6),  (7))  and  the  series  of  papers  by  ^ 
Agranovic  [(8),  (9),  (10)].  All  of  this  work  was  stimulated  by  the  pioneering  pape 
of  Kacenelenbaum  (11)  and  its  sequel  by  Voitovic,  Kacenelenbaum  and  Sivov  (12).  ■/  j 
There  also  is  an  important  paper  by  Agranovic  and  Golubeva  (13).  [A  word  of  cauticn 
must  be  injected  here.  In  a  relevant  paper  of  Golubeva  (14),  the  word  "proposition’ 
is  translated  as  "conjecture,"  which  one  must  admit  does  change  the  flavor.  Since! 
each  "conjecture"  is  followed  by  a  full  proof,  the  reader  must  be  alert  to"  the  ! 


translation  problems.] 


SINGULARITY  EXPANSION  METHOD (SEM) 


'/ 

) 

/  : 

/  i 

I 


The  Singularity  Expansion  Method  and  its  relationship  to  the  mathematical  theory  of 
scattering  has  recently  been  analyzed  by  Dolph  and  Cho  (15)  and  consequently  no 
exhaustive  detail  on  the  technique  is  required  here.  Instead,  a  summary  of  the 
results  of  that  work  will  be  given. 

‘  '  : 

The  Singularity  Expansion  Method  is  a  generalization  of  well-known  techniques  of 
linear  circuit  theory  in  which  the  singularities  of  a  transfer  function  are  used 
to  determine  the  transient  response  by  the  Heaviside  expansion  theorem.  In  electro 
magnetic  theory,  the  singularities  are  found  by  first  applying  a  two-sided  Laplace 
transform,  parameter  s,  to  the  Maxwell  equations  and  then  constructing  an  integral 
equation  for  the  scattered  field.  Complex  singularities  (sn)  appear  as  poles  of 
the  inverse  of  this  equation  and  are  determined  from  the  non-trivial  solutions  of 
the  corresponding  homogeneous  integral  equation. 


The  scattering  operator  is  a  unitary  operator  on  a  Hilbert  space.  Its  I.aplace 
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transform,  the  scattering  matrix,  is  analytic  in  the  right-hand  s-plane  including 
the  axis  and  is  meromorphic  in  the  left-half  plane.  It  has  poles  in  the  left-half, 
plane  at  those  values  of  s  for  which  there  exist  non-trivial  outgoing  solutions  of  ; 
the  reduced  Maxwell’s  equations  satisfying  the  boundary  conditions.  These  discrete 
values  of  s  are  complex  eigenvalues  of  the  exterior  scattering  problem.  Shenk  and  ; 
Thoe  (16)  have  established  a  one-to-one  correspondence  between  the  poles  of  the  j 
'  S-matrix  and  the  poles  of  the  integral  equation.  I  ‘  ! 

•  I  I 

Some  of  these  poles  occur  on  sheets  of  Riemann  surfaces  and  serious  difficulties  j 
arise  in  interpreting  their  meaning.  This  fact  can  be  deduced  from  the  study  of  ■ 
the  integral  equation.  Only  the  Green’s  Kernel  can  be  continued  analytically  into; 
such  "forbidden  domaxns".  fhe  resolvent  can  net.  j  j 

!  t  .  | 

While  it  is  difficult  to  relate  these  concepts  in  general  since  most  of  the  work  sn 

SEM  involves  a  formalism  in  which  neither ' spaces  nor  properties  of  the  .integral  "  j 
equation  are  given,  this  can  be  done  for  the  papers  of  Marin  (3).  He  uses  a  Hilbert 
space  consisting  of  tangential  currents  oh  the  surface  of  a  convex  body  to  discuss  ! 
the  magnetic  field  integral  equation.  He' deduces  a  Fredholm  integral  equation  of 
the  second  kind  from  this  and  uses  Carleman's  Fredholm  theory  for  the  determination 
of  the  natural  modes.  The  solution  of  the  corresponding  homogeneous  integral  1 

equation  contains  both  exterior  and  interior  resonances,  the  latter  being  purely  | 
imaginary.  .  i  |  •  i 

i 

To  relate  the  scattering  matrix  to  the  integral  equation  one  makes  use  of  its  ! 

representation  as  a  compact  Fredholm  integral  operator.  The  kernel  of  this  opera-'. 
tor  is  a  transmission  coefficient  arising  from  the  asymptotic  form  of  the  scattered- 
field.  Since  its  determination  involves  the  solution  of  the  scattering  problem,  1 
rather  than  use  the  above  representation,  it  is  simpler  to  use  the  one-to-one  ; 

correspondence  between  the  kernel  of  this  representation  and  the  solution  of  the  ! 
vector  wave  equation.  It  can  be  given  in  terms  of  a  vector  integral  equation  which 
is  more  general  than  the  magnetic  field  integral  equation.  Fredholm  theory  is  now’ 
used.  Recall  that  if  A  is  a  compact  operator  (such  as  the  above  mentioned  operator 
defining  the  integral  equation)  then  the  first  parf-  of  the  Fredholm  alternative 
states  that  A$  =  f  has  a  unique  solution  $  if  A<$>  =  0  has  only  the  trivial  solution. 
Here  the  first  part  of  the  Fredholm  alternative  yields  a  unique  solution  to  the 
general  equation  for  the  right-hand  s-plane  including  the  imaginary  axis.  The  ; 
analytic  Fredholm  theorem  for  compact  operators  which  is  given,  for  example,  in  j 

Reed  and  Simon  (17),  then  implies  the  analytic  and  meromorphic  properties  of  the  ; 

scattering  matrix  discussed  above.  ;  ,  /  i 

Vi 

•  .1  I 

For  the  magnetic  field  integral  equation, : one  must  use  the  second  part  , of  the  Fred¬ 
holm  alternative,  namely:  If  A<$>  =  0  has  non-trivial  solutions  then:  .  i 

,  •  |  > 

•  1  i 

(i)  =  0,  A*<^  .=  0,  where  A*  is  the  adjoint  of  A,  have  the'  same  finite  : 

number  of  solutions.  >  /  t  • 

'  i  /  I  ! 

(ii)  For  A<p  =  f  to  have  a  solution  <p^,  f  must  be  orthogonal  to  all  of  the  ! 

solutions  of  A*tp  =  0.  1  j 

J  i 

Since  the  solution  to  A<J>  =  0  can  be  added  to  then  clearly  the  process  is  not 
unique  and  no  one-to-one  map  between  the  scattering  matrix  and  the  integral  cquatio. 
exists  for  the  same  half-plane.  The  non-trivial  solutions  of  the  corresponding 
homogeneous  integral  equation  occurring  for  complex  {sn>  do,  however,  correspond  to 
poles  of  the  scattering  matrix  but  the  interior  resonances  corresponding  to  purely 
imaginary  {s  }  are  not  poles  of  the  scattering  matrix  nor  do  they  appear  in  the 
solution  fornthe  scattered  field.  They  are  consequently  method  dependent.  The 
exterior  resonance  corresponding  to  complex  {s^}  are  method  independent  and  intrin¬ 
sic  to  the  scattering  body. 


.  METHODS 

Since  a  one-to-one  correspondence  fails  to  exist,  it  is  conceivable  that  there  are  1 
complex  poles  of  the  scattering  matrix  which  are  not  given  by  the  solutions  of  the  ; 
homogeneous  magnetic  integral  equation  corresponding  to  complex  {sn}.  No  examples 
however,  in  which  this  occurs  are  known  and  thus  the  problem  remains  open.  | 

!  »  j 

Since  it  is  generally  believed  that  the  scattering  matrix  contains  all  observable  : 
information  about  the  scattering  process,,  the  above  observations  substantiate  the  ' 
claim  that  the  complex  singularities  are  intrinsic  to  the  scatterer  whereas  those  ; 

which  are  purely  imaginary  are  not.  i  * 

• 

Further  difficulty  arises  in  that  most  problems  treated  by  the  SEM  formalism  in- 
vole  /numerical  techniques  applied  to  matrix  equations  obtained  from  finite  approx¬ 
imations  to  Ljie  integral  equations.  Hence,  Kilters  -space  solutions  are  not  really, 
appropriate^  Fortunately,  the  results  can  be  given  in  terms  of  solutions  in  the 
Banach  space  of  continuous  functions.  In' this  space  the  Fredholm  theorem  as  given; 
by  Steinberg  (18)  can  be  used  to  discuss  the  magnetic  field  integral  equation  and  { 
to  establish  the  one-to-one  map  between  the  scattering  matrix  and  the  associated  ■ 
set  of  vector  integral  equations.  This  map  between  the  scattering  matrix  and  the  • 
integral  equation  has  been  carried  out  for  convex,  perfectly  conducting  bodies  only 
[see  Lax  and  Phillips  (19),  Brakkage  and  Werner  (20),  and  Dolph  and  Cho  (15)).  j 
Certain  other  aspects  of  the  SEM  formalism  have  also  been  investigated  by  Dolph  and 
Cho  (15)  and  serious  doubts  were  raised  regarding  those,  parts  of  it  where  integral 
equations  of  the  first  kind  are  used.  Also,  there  appears  to  be  a  confusion  in  the 
SEM  literature  between  the  eigenvalues  of  the  integral  equations  and  those  of^the  ’• 
vector  wave  equations  and  a  belief  persists  that  what  is  true  for  the  finite  matrix 
equations  holds  in  the  limit.  i  '  j 

Though  said  elsewhere,  it  seems  worth  repeating  that  areas  in  which  the  formalism  ’ 

needs  further  investigation  include:  :  ! 

» 

!  i 

I.  The  construction  of  variational  principles  useful  for  providing  estimates  for  ; 
the  location  of  the  poles  in  the  SEM  and  possibly  of  use  in  establishing  their  j 
existence  for  off-axis  poles  not  covered  by  the  known  Lax-Phillips  results  (21)  for 

the  scalar  case  and  their  generalization  by  Beale  (22)  for  the  electromagnetic  case 

•  /  ; 

»  •  *  J 

II.  An  investigation  of  integral  equations  of  the  first  kind  as  used  in  SEM  and  a' 

possible  justification  of  their  use  through  regularization  methods,  such  as  those  - 
of  Tihonov  (23).  *  '  j 


III.  The  extension  to  the  electromagnetic  case  of  theorems  sufficient  to  guarantee 

that  all  poles  are  simple.  This  should  include  generalizations  of  the  theorems  of 
Steinberg  (18)  and  Howland  (24) .  ;  ,  '■ 

'  f  f 

IV.  An  investigation  of  the  entire  functions  which  arise,  through  the  Mittag- 
Leffler  theorem, |  in  the  formalism.  This  should  include  the  determination  of  con-  - 
ditions  under  which  they  do  not  occur  and  their  explicit  form  when  they  do  occur. 


V.  The  creation  of  a  systematic  theory  of  the  asymptotic  contribution  of  branch 
lines.  As  can  be  seen  from  the  discussion  in  Dolph-Scott  (1),  such  a  theory  could 
have  implications  in  many  areas.  ; 

i  * 

t 

At  this  point  a  few  remarks  on  the  T-matrix  would  seem  to  be  in  order.  The  T-matri 
formalism  for  scattering  has  proved  to  be  an  efficient  way  of  obtaining  numerical 
results  in  a  number  of  complicated  problems  [see  Bolomcy  and  Wirgin  (25)).  It  was 
used  by  Waterman  (26)  in  1969  and  subsequent  publications  include  papers  by 
Peterson  and  Strom[(27),  (28),  (29)],  Strom  [(30),  (31),  (32)]  and  Varatharajulu 
and  Pao  (33).  However  there  is  a  relationship  between  the  T  and  S  matrices,  mmel> 
T  =  S  -  1,  as  stated  in  V.’u  and  Ohmura  (34),  and  physicists  in  the  past  have  used 
S  in  preference  to  T. 
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F.I GENMODE  EXPANSION  METHOD 


Baum  (2)  has  advanced  the  idea  of  synthesizing  transient  responses  by  means  of  the  : 
eigenmodes  of  integral  equations  of  the  first  hind  describing  the  system  response.  : 
Considerable  caution  must  be  exercised  in  this  approach  however,  since  as  Ramm  [(6) 
(7)]  has  pointed  out,  incorrect  results  can  sometimes  arise  due  to  the  non-self-  ; 
adjoint  property  of  the  operators  that  are  treated.  At  this  point  it  appears  j 

certain  that  the  theory  of  non-self-adjoint  operators  can  be  used  to  contribute  * 
substantially  to  the  understanding  and  limitation  of  the  formalism  of  EEM  (  and  ’ 
SEM) -  While  the  idea  of  using  non-self-adjoint  operator  theory  in  scattering  j 

problems  is  not  new - it  was  already  suggested  by  Dolph  (35)  in  1960 — it  does  not  •; 

appear  to  have  been  used  in  connection  with  SEM  and  EEM  in  the  English  literature. 
This  is  rather  surprising  since  this  theory  has  been  imployed  in  connection  with  | 
scalar  diffraction  problems  in  papers  translated  from  the  Russian  and  reprinted  in 
the  journal  Radio  Engineering  and  Electron  Physics  beginning  after  the 'formalism  : 
presented  in  the  paper  by  VortoviC,  Kacenelenbaum  and  Sovov  (12)  .  1  ■  i 

i  !  | 

While  this  latter  paper  contained  a  formalism  similar  to  EEM  for  both  the  scalar-  * 

and  electromagnetic  case,  including  dielectric  problems,  its  subsequent  interpretat 

at  least  in  translated  papers  known  to  the  writers,  in  terms  of  non-self-adjoint 
operator  theory  have  been  limited  to  the  scalar  problem.  In  the  simplest  case  one' 
attempts  to  construct  a  formal  solution  to  the  system  \ 


(A  +  k  )u  =  f 

u  =  0  ,  on  a  smooth  convex  scatterer  F, 


CD 


satisfying  the  radiation  condition  as  a  series  u  -  uq  +  EA^^x)  .  Here  uQ  is  the 
incident  field  and  (<Pn(x)}  are  the  eigenfunctions  of  a  compact~integrai  operator  ; 
A,  which  for  in  three  dimensions,  is  given  explicitly  by  1 


A$n  =  jGo(-~^n^dCT  * 


(x) 


(2) 

/ 


where  the  free  space  Green's  function  Gn  is 

^  I 
2 


J0  jx-y| 

I 


(3) 


l 


It  was  further  assumed  that  the  coefficients  in  the  above  expansion  could  be  determ 

by  the  Fourier  coefficient  formula  ,  /  1  l 

/  '  « 

/  :  ; 


Jr(u-W° 


m 


Ramm  (7)  interpreted  and  clarified  these  results  by  using  a  Hilbert  space  b  ( F) 
with  the  usual  Hermitian  inner  product 


(f.g) 


|f (x) g (x) do 

r 


(S) 
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Since  G„  is  real  symmetric  but  complex  valued,  this  implies  tliat 


(A$,t|>)  =  C <*> » Atp)  •  (6)  • 

» 

so  that  the  operator  A  is  non-self-adjoint.  !  ! 

Since  the  compact  operator  A  is  non-self-adjoint,  it  may  have  root  vectors  instead 
of  simple  eigenvectors.  That  is  for  a  given  X,  there  may  exist  an  integer  p  >  1  J 
such  that  (A-XI)^  =  0  for  some  4),  while  (A-XI)^  ^  0  for  all  q  <  p.  (In  the  matrix 
case  this  happens  when  non-simple  elementary  divisors  occur  and  requires  the  use  of 
the  Jcrd',n  nor*"51!  form  rather  than  the  diagonal  form  in  the  canonical  representatioi 
of  the  matrix.)  However,  Ramm  was  able  to  show  that  while  the  system  of  root  vector 
was  always  complete  in  I^CT),  the  simple  form  of  the  coefficients  given  above  would 
only  occur  if  the  surface  r  were  such  that  A  defined  over  it  was  a  normal  operator, 
that  is  AA*  =  A*  A  (this  condition  is  necessary  and  sufficient  for  A  to  be  a  diagonal 
operator)'  The  normalit>  of  the  operator  can  be  tested  as  follows:  For  example,  i 
consider  the  operator  A:  H->H,  where  J  [ 


A  is  normal  provided  the  integral 


vanishes,  a  condition  that  Ramm  was  able  to  show  was  true  if  F  is  a  sphere.  In 
complicated  problems  it  may  have  to  be  tested  by  direct  computation.  ! 

I 

* 

If  the  operator  A  is  not  normal,  considerable  complexity  arises  in  the  theory  and 
computational  schemes.  The  problem  is  that  in  general  one  has  to  work  in  an  in¬ 


finite  dimensional  analog  of  the  Jordan  normal  form,  since  A  fails 'to  be  dia-  ; 

gonalizable.  The  lemma  of  Schur  [see  Gohberg  and  Krein  (36)]  [states  that  any 
completely  continous  operator  A  mapping  a  Hilbert  H  space  onto  itself  can  be  repre¬ 
sented  in  a  triangular  form.  Specifically,  there  exists  an  orthogonal  basis  to .  of. 
H  such  that  :  I  .  J  i 


Ato.  =  X 

l  .  .  Ji  i 

J  i<J 


where  a..  -  (Ato., to.) 
of  A.  31  J  J 


.  •  » 

X.,  (3)  [denoting  inner  product  and  X.  being  an  eigenvalue 

3  i  J 


There  are  many  areas  which  need  detailed  mathematical  investigation.  For  example 
VoTtovic,  Kacenelenbaum  and  Sivov  consider  several  variational  principles  which 
produce  stationary  solutions.  These  are  reminisiccnt  of  those  of  Schwinger  and 
McFarlane,  the  latter  occuring  in  the  problem  of  anomalous  propagation  through 
the  atmosphere.  Attempts  to  make  mathematical  sense  of  these  occur  in  Dolph  (37), 
Dolph  and  Ritt  (38)  and  Dolph,  McLaughlin  and  Marx  (39).  The  max-min  characterizat 
of  the  last  paper  unfortunately  depends  upon  the  dimension  of  the  approximating 
space  and  so  it  badly  needs  reformulation.  The  only  other  pertinent  work  it  seems 
to  be  that  of  Morawetz  (40)  where  as  discussed  in  Dolph -Scott  (1)  a  variational 
principle  is  given  implicitly. 
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In  connection  with  variational  principles,  the  following  facts  should  be  mentioned. 
The  variational  approach  given  in  (5)  has  been  briefly  reviewed  by  Doiph  in  (41).' 
Voitovic,  Kacenelenbaum  and  Sivov  introduced  the  e-approach  by  considering  an 
approximate  scattering  problem  (the  e-problem)  and  the  corresponding  systems  of 
equations  1 

I  i 


.2 

k  eu  = 

2  1 
K  u  =i 


-  u"!< 


in  V 


leading  to  a  miminization  of  the  functional 


(Vu)  2dV 


2 

udV 


A  rigorous  justification  of  such  techniques  still  remains  an  open  problem.  Similar 
variational  principles  were  given  by  J.  Schwinger  and  H.  Levine  (42)  .  [Also  see  ) 
Kate  (43),  and  Doiph  (37)].  As  in  the  case  discussed  by  Doiph  and  Ritt  (38),  or 
Doiph  (37j,  it  cm  be  conjectured  that  the  real  and  imaginary  parts  of  the  unknown 
function  near  the  stationary  point  lie  on  a  saddle-like  surface,  but  the  orientation 
of  such  a  saddle  is  unknown.  '  •  ! 
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ZEROS  OF  THE  INPUT  IMPEDANCE  FUNCTION  FOR  THTN 
CYLINDRICAL  AND  PROLATE  SPHEROIDAL  ANTENNAS 

C.  T.  TAI 

1.  Introduction 

The  input  impedance  function  of  the  thin  cylindrical  antenna 
has  been  studied  by  Hallen  [1],  Schelkunoff  [2],  King-Middleton  [3], 

Tai  [4]  and  others.  Although  the  methods  used  by  these  authors  are 
different  from  each  other,  yet  they  yield  the  numerical  results  for 
the  input  impedance  in  the  (real-valued)  frequency  domain  which 
agree  relatively  closely  not  only  with  each  other  but  also  with  ex¬ 
perimental  results. 

The  main  objective  of  our  study  is  to  numerically  compute  the 
distribution  patterns  of  zeros  of  the  input  impedance  functions  in 
the  complex  frequency  domain  of  thin  cylindrical  and  prolate  spher¬ 
oidal  antennas,  and  study  their  implications.  This  is  carried  out 
by  the  use  of  Schelkunoff ' s  perturbation  method  [2]  with  a  slight 
modification.  The  Schelkunoff ' s  method  is  essentially  a  heuristic 
one  and  so  is  our  present  study  because  we  adopt  the  method  formally 
without  establishing  a  rigorous  mathematical  justification. 

In  section  2  the  input  impedance  functions  are  derived  for  a 
cylindrical  antenna  and  a  prolate  spheroidal  antenna  and,  in  section 
3,  the  numerical  results  are  presented  for  zeros  of  these  functions 
in  the  complex  frequency  plane.  In  section  4  we-  offer  a  concluding 
remark  on  our  numerical  results. 

2.  Input  Impedance 

In  this  section  we  derive  the  input  impedance  functions  for  thin 
cylindrical  and  prolate  spheroidal  antennas  based  on  the  Schelkunoff ' s 
method.  For  convenience  of  the  reader,  we  will  show  the  essence  of 
the  method  in  some  details  for  the  cylindrical  case;  for  the  prolate 
spheroidal  case,  we  will  merely  show  the  result  since  the  same  procedure 
is  used. 

2.1  Cylindrical  antenna 

Consider  a  cylinder  of  lengtn  21  and  radius  a.  By  regarding  it 
as  a  nonuniform  transmission  line  with  the  source  located  at  the 
center  of  the  cylinder,  we  define  the  line  impedance,  line  admittance, 
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and  characteristic  impedance  at  r,  0  <■  r  <-  l,  as  following: 

no  2r 

Z(r)  =  jk  ^  In  i~~) 

Y  (r)  =  jk-rr/n  In  (--)  , 

U  d 

and 

no  2r 

K(r)  =  -2  In(^) 

7T  a 

where  n  =  120tt  ohms  and  k  =  w/u  C  T 
o  o  o 

We  then  define  the  average  characteristic  impedance  of  the  non- 
uniform  line  as 

Ko  “  ir  I  /*  a)dr'  (1) 

o 

where  a  is  some  positive  real  constant  number  yet  to  be  determined. 
Carrying  out  the  integral,  we  get 


where  3 
eristic 


v  =  °  o  n  ( 2^.  1\ 

Ko  tt  £n(a  F  ' 


(2) 


=  e/a.  If  we  set  3  =  e,  (2)  reduces  to  the  average  charact- 


impedance  defined  by  Schelkunoff: 


K  1 

°J3=e 


"o  .  ,21, 

IT  ln  <«>  ■ 


which  can  be  written  in  terms  of  the  bicone  parameter  0Q  =  a/I  [5], 
the  half-angle  of  the  bicone,  as 


k  1  =  ^ 

°J  o  tt 

J  0=e 


*n(le> 


This  shows  that  the  average  characteristic  impedance  of  the  cylinder 
corresponds  to  the  characteristic  impedance  of  a  bicone  whose  half¬ 
angle  is  e0Q  -  2.7  0  .  On  the  physical  ground,  it  seems  more  appro¬ 
priate  to  have  the  bicone  inscribed  in  cylinder.  It  is  for  this 
reason  that  we  introduced  a  new  averaging  procedure  in  the  form  of 
(1).  If  3  =  1  is  chosen,  (2)  becomes 
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K  =  —  Z  n(~)  =  Z  (0  ) 
O  IT  0  '  C  O 


2c(eo)  bein<3  the  characteristic  impedance  of  a  biconical  antenna 
with  the  half -angle  dQ  =  a/Z. 

The  average  line  impedance,  Zq,  can  also  be  defined  similarly. 

Thus 


H  rZ  ~ 

=  -¥  I  J  *n(ir  a)dr 


jkK  . 
J  o 


Following  Scheikunoff,  we  set 


(4) 


and 


k 


/Z  Y  . 


o  o 


(5) 


Note  that  Yo  is  not  che  average  admittance  in  the  sense  of  (4).  We 
now  denote  the  measure  of  nonuniformity  of  the  line  by 


2(r) 

=  Z  (r)  -  z0. 

(6) 

O 

A 

A 

9(  r) 

=  Y(r)  -  Yq. 

(7) 

Then,  we  can  write 


z(r)  =  jk  ^  Zn(j)  ,  0  £  r  <  Z.  (8) 

This  "perturbation"  of  line  impedance  is  awkward  because  it  tends 
to  -  j=°  as  r  ->  0.  This  is  a  consequence  of  the  definition  of  the 
line  impedance  as  given  previously,  and  is  one  of  the  mathematical 
difficulties  in  the  Schelkunof f 1 s  method.  We  then  assume  that  the 
nonuniform  line  varies  in  such  a  way  that  the  following  restriction 
holds  throughout  the  line: 


&l£>  =  _ 

2  (r )  Zo 


<9) 
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This  condition  is  necessary  for  linearization  of  the  theory. 

By  the  nonuniform  transmission  line  theory,  the  line  voltage 
and  current  at  r,  0  <  r  <.  I ,  are  given  by 


and 


where 


V  {r )  -  VQ(r) 


-k  l 


+KqI (t) 


r 

2(t)  [  jV  (t)  sink  (r-t) 
cosk  (r-t)  ]dt 


(10) 


I  (r) 


Vo(r) 


q (r) 


+V(t) 


=  V  coskr  - 
o 


=  I  coskr  - 
o 


2 (t) [ jKQI (t) sink (r-t) 


cosk  (r-t)  ]  dt , 

(ID 

jK  I  sinkr, 
o  o„ 

(12) 

V 

j  sinkr, 

(13) 

o 


IQ  =  I  (o)  ,  VQ  =  V(o)  . 


VQ(r)  and  IQ(r)  denote,  respectively,  the  line  voltage  and  current 
of  a  uniform  transmission  line.  Eqs.  (10)  and  (11)  are  coupled 
Volterra  integral  equations  of  the  second  kind  for  I (r)  and  V(r). 

Since  such  Volterra  integral  equations  in  possess  unique  solutions 
and  the  operators  are  contractive,  one  can  apply  Piccard's  iteration 
to  obtain  a  convergent  solution.  Thus,  we  formally  write  the  solutions 
of  (10)  and  (11)  as 


00 


V  (r ) 

=  £  V(r) 

n=0 

oo 

(14) 

I  (r) 

=  E„  In(r) 
n=0 

(15) 

For  n  =  0,  VQ(r)  and  IQ(r)  are  given  by  (12)  and  (13);  for  nil. 


Vn(D 


—  f  z  (t )  [  j  V  .  (t)  sink  (r-t) 
Ko  Jo  n_1 


I 

n 


(r) 


+  KQIn_^(t)  cos  k  (r-t)Jdt,  (16) 


[ jK  I  .  (t ) sink  (r-t ) 
o  n- 1 


+  V 

n- 


1 


(t ) cosk (r-t) ] dt 


(17) 


Unfortunately,  it  becomes  impractical  to  compute  the  iterates  be¬ 
yond  the  first  and  this  is  a  major  difficulty  of  this  theory.  We 
are  then  forced  to  accept  the  "approximate"  solutions  as 


V  (r )  =  Vq  (r)  +  Vx  (r)  , 


I  (r)  =  IQ  (r)  +  I1 (r) . 


Following  Schelkunoff,  we  set  the  terminal  impedance  at  r 


Z(k£)  = 


V  (k£ )  .  Vo(k^  +  Vl(k£) 


I  (k£)  IQ  (k£)  +  I  (k£) 


(18) 

(19) 


l  as 
(20) 


We  note  that  V^(k£)  and  (k£)  in  (20)  are  the  first  iterated  solu¬ 
tions  of  the  Fredholm  integral  equations  with  r  replaced  by  l  in  (10) 
and  (11).  As  such,  unlike  the  Volterra  integral  equations,  we  no 
longer  have  a  priori  assurance  of  the  uniqueness  of  the  solutions, 
nor  of  the  convergence  of  the  iterated  solutions.  Even  if  we  assume 
that  tne  iterated  solutions  of  the  forms  (18)  and  (19'  are  acceptable 
approximate  solutions  for  k  >  0,  (i.e.,  for  the  real-valued  frequen¬ 

cies),  i+-  is  possible  that  such  an  assumption  may  not  be  valid  when 
analytically  continued  in  the  complex  frequency  plane  as  we  must  even¬ 
tually  do  for  the  computation  of  zeros  of  the  input  impedance  function. 
In  spite  of  these  unsettled  mathematical  questions,  we  formally  pro¬ 
ceed  and  postulate 

K  2 

z(kt>  -  rW  ’  1211 

where  a 

(k)  =  K  2Y  (k£)  .  (22) 

cl  O 


Y(k£)  [5]  represents  the  terminal  admittance  of  a  biconical  antenna 

with  which  our  cylindrical  antenna  is  associated.  The  configuration 
of  the  biconical  antenna  depends  on  the  choice  of  the  value  of  8- 
Evidently,  the  postulations  (21)  nd  (22'  are  heuristic.  Equating 

(20)  and  (21) ,  we  can  find  the  i  impedance  function  of  the  cyl¬ 

indrical  antenna  at  r  ~  0.  This  is  the  Schelkunof f ’ s  theory  of  the 
cylindrical  antenna.  Actually,  Schelkunoff  approximated  (20)  as 
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Z(k  i)  = 


Vo(k<l)  +  Vx(k£) 

-^kT) 


At  this  point,  we  depart  from  the  Schelkunof f 1 s  theory  and 
define  a  complex  variable  s  by 


s  =  jk5. 


(23) 


and  formally  replace  k£  in  (20)  and  (21)  by  -  js.  After  some  lengthy 
manipulation,  we  obtain  the  desired  complex-valued  input  impedance 
of  the  cylindrical  antenna  in  s  domain  and  write 


where 


e’c  ( s )  =  f  ( s ; 


(c)  Fc(s) 

z:  is)  =  k  -  -  ■  , 

m  o  G  (s) 


(24) 


1  -r-  [eS  E  ( 2  s )  +  e  SE(-2s)] 


20  -  4i,n3 


-  ( e S E  (2s)  -  e  sE(-2s)  )  , 


2  SI 


(25) 


Cc(s) 


g  (s) 


20  -  4S,n3 


4£n$sinh (s) 


[ esE ( 2s )  -  e  SE(-2s)] 


_  y  (s)  ( 
20  ) 


4£n0cosh  s) 


( 


[eSE  (2s)  +  e  SE  (-2s) ]>.  (26) 


f  (s)  =  J-  [1  +  y  (s)  ]  [1  +  f(s)e  ZS]  , 


g(s)  =  [1  i-  y  (s)  ]  [  1  -  f(s)e  2s]  , 


'I' 


(27) 

(30) 


y  (s) 


2s 


-2s. 


20  l  in  (2)  •  [e"°  -  e  ‘"'3 

(+  (2  +  e2s)E(2s) 

+  e"2SE(-2s) 

-  e25E(4s) 


(31) 
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and 


E(±s) 


1 


2u 

-  e 


u 


du 


SI  =  2£n(|-)  . 

°o 


(32) 

(33) 


The  input  impedance  of  the  biconical  antenna  of  the  half-angle  0^ 
is  given  by  (cf[5j) 


z.(b)  (s: 

in 


Z  (0  ) 
c  o' 


f  (s  ) 
g  (s) 


z  (0  ) 

c  o 


1  +  r (s)  e 
i  -  r (s) 


-2s 


-2s 


(34) 


The  extent  of  the  validity  of  the  input  impedance  function  (24) 
can  be  tested  by  the  inicial-  and  final-value  theorems.  Our  analysis 
shows  that,  with  the  applied  unit  step  voltage,  the  final-value 
theorem  is  satisfied  but  not  the  initial -value  theorem.  In  fact, 
we  get 


£im  sK  I  (s )  •>  -  1 , 

j  s  |  ->•  oa 

which  is  clearly  absurd.  This  finding  cast  some  serious  doubt  about 
the  numerical  values  of  the  zeros  of  the  input  impedance  function.* 


2.2  Prolate  spheroidal  antenna 

We  now  consider  a  thin  prolate  spheroidal  antenna  inscribed  in 
the  cylinder  described  in  section  2.1.  By  the  same  procedure,  we  can 
obtain  the  formal  expression  for  the  input  impedance  function  in  s. 
The  result  is 


'./here 


.(s)  (s)  =  K 

‘in  (s)  Ko  G  (s) 
s 


(s) 


Fc(s)  + 


D, (s)  ,  , 

_i .  y (s)  •  p  ,  i 

n  -  2 Znfr  +  SI  jEl  (  } 


(35) 

(36) 


*The  approximate  impedance  function  given  by  eg.  (8)  by  Hallen  (1) 
also  fails  to  meet  the  initial-value  theorem. 
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Gg  (s)  -  Gc  (s) 


jE1 (s) 
ft  -  2Jln8 


y  (s) 
ft 


Dj_  (s)  , 


D^(s)  =  -  sinh  (s)  sinh  ( 2s )  •  Jin  (2 ) 

-  i  cosh (s)  [e2sE (2s)  +  e“2sE(-2s)] 
+  i  [e3sE (4s)  +  e“3sE(-4s)] , 


jE^(s)  =  -  cosh  { s )  sinh  (2s)  •  Jin  (2) 

-  j  rinh(s) [e2sE(2s)  +  e~2sE(-2s)] 
+  J  [e3sE(4s)  -  e"3sE(-4s)  ]  . 


(37) 


(38) 


(39) 


3.  Numerical  results 

We  have  numerically  computed  zeros  of  t  he  complex-valued  input 
impedance  functions  (24)  and  (35)  in  the  s-plane  with  8=1  for 
a/Z  =  0.01  and  0.1,  respectively  and  the  results  are  shown  in  Figs. 

1  and  2.  In  Fig.  1,  we  also  show  corresponding  zeros  of  the  bicon- 
ical  antenna  for  a  direct  comparison.  As  seen  from  this  graph,  the 
zeros  of  the  input  impedance  function  for  the  cylindrical  antenna  are 
only  slightly  displaced  from  these  for  the  inscribed  biconical  antenna. 
For  the  purpose  of  attempting  to  understand  the  implications  of  the 
numerical  results  of  our  present  study,  we  also  include  in  the  graph 
complex-valued  poles  of  the  axial  current  density  function  on  the  cyl¬ 
inder  of  the  dimension  a/Z  -  0.01.  The  data  shown  in  Fig.  1  are  ob¬ 
tained  by  Tescne  [6]  by  means  of  the  singularity  expansion  method  or 
the  integral  equation  method;  Fig.  1-1  shows  the  distribution  pattern 
of  these  complex  poles.  The  SEM  complex  poles  correspond  to  complex¬ 
valued  zeros  of  the  input  impedance  function.  We  see  from  Fig.  1  that 
the  distrioution  pattern  of  our  results  and  that  of  the  SEM  are  en¬ 
tirely  different.  In  particular,  the  zeros  of  the  input  impedance  based 
on  the  Schelkunof  f '  s  theory  are  (located  only  on  two  "layers",  whereas 
the  Tescne' s  poles  occur  in  mul fi-layers  of  similar  trajectory.  it 
should  be  emphasized,  however,  that  the  nature  of  our  present  problem 
and  that  of  the  SEM  are  dicferent:  our  cylinder  is  a  center-fed  trans- 
mi  cting  antenna,  whereas  the  cylinder  in  the  SEM  analysis  is  a  scatter- 
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ing  body.  Therefore,  there  is  no  basis  to  expect  two  sets  of  data 
to  have  similar  variations. 

In  Fig.  2,  we  present  zeros  of  the  input  impedance  function  of 
the  prolate  spheroidal  antenna  of  the  dimension  a/l  =0.1  and  the 
corresponding  SEM  results  obtained  by  Marin [7];  the  distribution  pat¬ 
tern  of  the  latter  is  shown  in  Fig.  2-1.  The  zeros  of  the  input 
impedance  function  of  our  prolate  spheroidal  antenna  occur  again 
only  on  two  layers,  but  the  layers  now  cross  each  other.  Such  a 
crossing  of  two  layers  also  occur  in  the  case  of  a  biconical  antenna 
when  the  half  angle  0q  is  relatively  large.  Again  our  numerical 
results  and  those  of  the  SEM  are  widely  different. 

1.  Conclusion 

The  objective  of  our  present  study  was  to  shed  some  light  on 
the  transient  behaviors  of  thin  cylindrical  and  prolate  spheroidal 
antennas  by  the  use  of  the  classical  pole-zero  method  of  the  transfer 
function.  Unfortunately,  however,  we  do  not  feel  that  we  can  draw 
definite  conclusions  about  the  implication  of  the  numerical  results 
for  the  zeros  of  the  input  impedances  of  the  antennas,  mainly  because 
v/e  are  not  certain,  as  discussed  in  section  2.1,  whether  or  not  the 
formal  extension  of  Schelkunof f ' s  theory  to  the  complex  frequency 
domain  is  meaningful,  even  though  Schelkunof f ' s  theory  is  generally 
believed  to  yield  the  input  impedance  function  of  a  thin  antenna 
structure  in  the  sense  that  the  numeral  results  agree  reasonably 
well  with  experimental  ones  in  the  real  frequency  domain.  Finally, 
we  remark  that  two  similar  sets  of  zeros  of  the  input  impedance  func¬ 
tion  in  the  complex  frequency  domain  such  as  those  for  the  cylindrical 
and  biconical  antennas  in  Fig.  1  do  not  necessarily  imply  that  the 
impedance  functions  in  the  real  frequency  domain  are  close  to  each 
other.  For  example,  the  maxima  and  the  minima-  of  the  real  and  imag¬ 
inary  parts  of  the  input  impedance  of  the  biconical  antenna  and  those 
of  the  cylindrical  antenna  in  Fig.  1  are  approximately  2  to  1  in  ratio. 
Conversely,  two  input  impedance  functions  which  agree  arbitrarily  closel 
with  each  other  in  the  real  frequency  domain  may  have  quite  different 
sets  of  zeros  in  the  complex  frequency  domain. 
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Fig.  l:  Zeros  of  input  impedance  function 
cylinderical  antenna  for  a/l  =.01 
based  on  the  perturbation  method. 
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Zeros  of  input  impedance 
function  of  a  prolate  spheroidal 
antenna  for  a/Z  =  0.1  (0  =  5.73°) 

based  on  the  perturbation  method. 


